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Two simple proofs are presented for the first order virial expansion of the self-energy of a particle 
moving through a medium, characterised by temperature and/or chemical potential(s) . One is based 
on the virial expansion of the self-energy operator itself, while the other is on the analysis of its 
Feynman diagrams in configuration space. 
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I. INTRODUCTION 

More than a decade ago, Leutwyler and Smilga Q considered the problem of mass-shift and damping rate of a 
nucleon propagating through a heat bath. As one of the methods to deal with the problem, they wrote (the first term 
of) the virial expansion of the nucleon self-energy, which relates it to the pion-nucleon scattering amplitude in the 
forward direction. At low pion density the first order formula yields the dominating contribution. The advantage of 
such a formula is that one can use the experimental data to compute the shift in self-energy without relying on any 
theoretical structure. Further the formula is simple enough to suggest a generalisation to other hadrons in different 
media. Accordingly a number of authors have used these relations to explore the properties of hadrons in such media 

tfl- 

In the context of quantum electrodynamics, Barton pi took a physical approach to derive such formulae, following 
an earlier suggestion by Feynman j(3| . The refractive index of a medium for a beam of particles of a definite frequency 
differs from unity due to the scattering of the incident beam by particles of the medium. The change in frequency 
(energy) of the particles may be related to this change in refractive index. In this way Feynman related the self- 
■ energy of the electron to its forward scattering amplitude with the virtual quanta in vacuum. Barton extended these 
' considerations to the black body radiation. 

Recently the complete expansion of the self-energy in the medium in powers of the distribution function has been 
proven in perturbation theory [Q. In effect, it is a reordering of the contributions in perurbation theory, which is an 
expansion in powers of coupling contant, to one in powers of the distribution function. The analysis is carried out in 
imaginary time formulation of field theory in the medium. 

Here we present two simple proofs of the first term in the expansion of the self-energy. One is a generalisation of 
the method used in Ref. |f| and is based directly on the virial expansion of the ensemble average of any operator. The 
other is in perturbation theory, based on Feynman diagrams in configuration space. 

We shall work here in the real time formulation of field theory at finite temperature and density. Sees. II and III 
describe the two proofs of the virial expansion. Sec. IV discusses the relevance of such formulae and possibility of 
extending the proofs to higher order. 



II. VIRIAL EXPANSION 



Following Weinberg ||, we introduce a compact notation for different types particles and their fields. Let a(p, a, f) 
and a' (p, a, f) be the destruction and creation operators for a particle of species / having momentum p and spin 
projection a; then for example, 



\p,a,f) =<J(p,a,f)\0). 



They satisfy the commutation (anticommutation) relations (u p = \/ p 2 + m 2 ), 

[a(p,a,/),at(pV,/'] T = (2tt) 3 2u p S 3 (p - p')5^S fr . 

The field operators are denoted by ipi(x), where the index I denotes not only the field type but also runs over its 
components. Then ipi(x) may be expanded as 
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/ c denoting the antiparticle of the species /. The coefficient functions ui and vi depend on the spin of the particle. 
We have in mind the three types of fields, namely, the scalar field 4>{x) for which u; = u; — 1; the Dirac spinor field 
ip(x), for which they are the normalised Dirac spinors, u(p, a)u{p, a) = —v(p, a)v(p, a) = 2m and the vector field A^, 
for which they are the polarisation vectors, t^p, a), e*(p, <r)e M (p, a') = —S arT i. Below we shall in most places suppress 
the variable a of ui, vi, a and . If there is an integration over a 3-momentum, a summation over the corresponding 
a as in Eq (2.1) will also be implied. 

In the real time formulation of field theory in medium, all Greens functions, in particular the self-energy, acquire a 
2 (g> 2 matrix structure. Since we confine here only to contributions linear in the distribution function, it will, however, 
suffice to consider only the 11-component of these matrices. Below we shall drop the 11 index []. 

Let us begin with the self-energy (p) of a particle of type / and momentum p in vacuum. For our purpose, we 
write this amputated, two-point Green's function as 

- i{2k)H\p' - p)uifa m ( p{p) Ul {p, /) = {0\a(p>, f)(S - l)at(p, /)|0), (2.2) 
where the .S-operator is given by the familiar time ordered expression, 

Q rp^L J d^xCj (x) 

Cj(x) being any interaction Lagrangian built out of the fields ipi(x). Because S?^(p) is one-particle irreducible by 
definition, it is understood here and below that we only retain such diagrams in matrix elements like the one on the 
right hand side of Eq.(2.2). The corresponding self-energy in the medium will be denoted by Hf{p,u), where u M is 
the 4- velocity of the medium Actually we shall work in the rest frame of the medium (u° = 1, u = 0). Then we 
may write 

- i (27r) 4 <5 4 (p'-p)u / (p,/)S / (p, U K(p,/) = (a(p'J)(S-l)rf(p, /)), (2.3) 

where (• ■ •) denotes ensemble average: For any operator O, 

(O) = Tr(pO), p = e- fj{H -^ } /Tre-^ B ~^\ (2.4) 

Here H is the Hamiltonian, (5~ x the temperature. For illustration we include a chemical potential /i for a fermionic 
species with the corresponding number operator Af. 
The ensemble average (O) admits a virial expansion 

(O) = (0\O\0)+J2 J {27T)32ujk n r (u, k )(k,f\0\kJ>) + ---, (2.5) 

where the sum over /' runs in general over the species of particles present in the medium. The dots represent terms 
of higher order in the distribution function n. The latter are given by the familiar expressions, nsi^k) = (eP u}k — 1) _1 
for bosons and rip(uik) — {e^^ kTIJ ^ + 1) _1 for fermions and antifermions respectively. Applying Eq.(2.5) to the left 
hand side of Eq.(2.3), we get for £/(p, u) = Y,f(p,u) — Sv'(|>), 

(2 * ^ /'KXfc, S'W, f)(S - l)a){p, f)\k, /')■ (2.6) 



1 The self-energy function E which actually shifts the pole position of the propagator is related to the components in a 
simple way E3|. With En these relations are, ReT,(p) = ReT,n(p) and JmE(p) = almEn(p), where a = tanh{f3po/2) for 
bosons and a — coth{/3(po — A t )/2} for fermions (p > 0). Thus En deviates from E by exponential corrections at small 
temperatures. Further the factor a actually cancels on rearranging the distribution functions within the integral, at least for 
1-loop contributions [p"l| . 

2 No confusion should arise from the use of same u in both u M and Ui(p, /). 
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The matrix element on the right in this equation is recognised to be the amplitude T for scattering of a /-particle of 
momentum p with a /'-particle of momentum k in the forward direction, 

(*, f'\a(f, f)(S - l) a t(a f)\k, /') = i (2^) 4 (S 4 (p' -0)T ff ,(p, k). (2.7) 
We thus get the virial expansion for the self energy to first order, 

/d 3 k 
J^3^ n f'( UJ k)T fr (p,k), (2.8) 

where an average over the polarizations of /-species and a sum over the polarizations of the different /'-species are 
understood. 

We now apply this formula to two cases of interest. The first one is that of a nucleon in a medium. Its complete 
propagator in the medium is m — £/). For the nucleon at rest (p = 0) in the rest frame of the medium, we can 

write £/ = a + 67 , getting w(p)£ fu(p) — 2m^(a + b). Also the complete propagator simplifies to 



Pa - (mN + a + b) 2 
Thus the new pole position is given by fl, 



i 1 v-^ f d 3 k 

fh N --r N = mN - — Y,J -^^n{^)T Nr [p,k). M 



/ 

The other case we consider is that of a vector meson in the medium. Let us rewrite Eq.(2.8) as 

d 3 k 

^ / 72 

f 

The free propagator for the vector meson is 



1 f rflk 

-]>>;>, a)n^, U )e,(p, a) = £ J ^—n f ,(u; k )T vf ,(p,k). 



(2.10) 



D$(P) = ~9,u + ^-L^ (2.11) 

To sum the series of 1-particle reducible insertions of the polarisation tensor, we have to decompose the latter in terms 
of kinematic covariants, 

%u = P^t + Q^u (2.12) 

which we choose as 

2 4 

p _ 1 VjlVv _ P_ ~ ~ p. P_ ~ ~ ;„ ,ql 

p z p z pr 



where co = u ■ p, p = y 1 ui 2 — p 2 and = — ujp^jp 2 . The covariants are free from singularities at p 2 = 0, but at 
p = 0, there is a constraint on the two amplitudes, 

n t ( P o,p = 0)=p 2 IL l ( Po ,p = Q). (2.14) 

Then the full propagator becomes 



u »» ~ P 2 _ m * v _ n t V - m2/(1 + n/ ) m 2 v 



Thus in general the transverse and the longitudinal components suffer different shifts in the position of the pole. But 
for the vector meson at rest (p — 0) the two shifts coincide, because of the constraint equation. We then get the same 
formula for the pole shift as Eq.(2.9) with the subscript N replaced by V. 
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III. PERTURBATION EXPANSION 



We now attempt a perturbative proof of Eq.(2.8). We begin by expanding the S-operator in Eq.(2.3) in the familiar 
perturbation series, 

-i(2w) 4 5 i (p'-p)u l (p,f)Z f (p,u)u l (p,f) = J2 jvT / d i x 1 ---d*x N FW(x 1 ---x N ), (3.1) 

N=l ' J 

where 

FW( Xl ... Xn ) = ( a (p', f)T{C x {xi) ■ ■ ■ d{x N )}a\p, /)} paircd . (3.2) 

The subscript 'paired' indicates that is the sum of all connected terms obtained by pairing (contracting) all the 

operators in it in all possible ways. In other words, it represents the sum of all connected Feynman diagrams in the 
iVth order. In the following we indicate a pairing by thick dots as superscript. The pairing between a creation or an 
annihilation operator of a particle with a field operator are given by, 

a'(p'JW(x)=e l P'-*u l (p'J) 

^(x)a^'(pj) = e^u l (pj). (3.3) 

We may choose to work out the self energy of a particle. But the medium may contain antiparticles. So we also note 
the contractions, 

a'(p'J c W(x)=e lp '- x vi(p'J) 

tiixtf'fa n = e- ip - x m /)• (3.4) 

Finally the pairing of two field operators results in the free propagator in the medium p2| , 

tffrXiv) = (TMx)4> m (y))- (3.5) 

A free propagator in the medium differs from that in vacuum if there are like-particles in the medium. In that case, 
it has an extra term containing the distribution function n of the particles in the medium and a mass-shell <5-function. 
Thus, isolating this term amounts to putting the internal line on mass-shell, i.e. opening the propagator into two 
external lines. When the ^-function is integrated out, this n-dependent piece in the propagator, to be denoted by a 
subscript n, becomes, 

(TM^m(y))n= J {2 f )3 k 2cJk {±nj(cj k )e- ik < x -^ Ul (k,f)u m (k,f) +nj(w k )e ik < x -^vi(k,f)v m (k,f)] , (3.6) 

where + (-) sign before the first term holds if the species / is a boson (fermion). The distribution functions n + and 
rT coincide if there is no chemical potential. 

Let us assume, for simplicity, that only one species of particles has its free propagator altered in the medium by 
the additional term (3.6). We single out this field from the compact notation ipi(x) and call it xi( x ) and denote its 
species by /' 0. We now wish to collect all the linear contributions in n by opening each of the x-propagators in turn 
in each of the Feynman diagrams for the self energy in the N th order of its perturbation expansion (Fig. 1). To 

this end we consider the sum of Feynman diagrams, to be denoted by F^\x\, ■ ■ ■ ,xn), containing a x-propagator 
between any two vertices, say at Xi and Xj, 

Fjf\xi, ■■■,x N ) = {a{p' f)T{d{ Xl ) ■ ■ ■ Ci'{xi) ■ ■ ■ Ci m (xj) ■ ■ ■ d(x N )}a\p, /)) paired , (3.7) 

where we explicitly indicate the pairing of the two fields at Xj and Xj. 

Consider first the indicated pairing in Eq.(3.7) before any other pairings are carried out, ie, keeping all other 
operators in their respective positions in the T-product. By a certain number of interchanges of the field operators 



3 For self-conjugate species like ty°, the two terms in Eq.(3.6) contribute to the same amplitude, the direct and crossed 
processes being identical. 
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we bring them together to form (xi( x i)Xm{ x j)}- Then we extract its n-dependent contribution given by Eq.(3.6) ^. 
Since we restrict to linear contributions in n, we set all other in-medium propagators to their vacuum values. 

In order to interpret the resulting expression in terms of a two-body scattering amplitude, we use (3.3) to rewrite 
the coefficient functions for /' species in {xi( x i)Xm( x j)) as pairings, 

e- ik < x *-^ Ul (k, f)u m (k, /') = x'{xi)<J\k, f)a"(k, f)x"( x j)- (3-8) 
For the antiparticle species f' c , we use Eq.(3.4) to write an analogous equation, 

e ik < x <-^ Vl (k, f)v m (k, /') = a'(k, f c )xK x i)x"( x iW"(k, f' c )- (3.9) 

Now the x fields can be brought back to their original positions forming again the complete vertices Ci{xi) and Ci(xj). 
We also bring the creation and the annihilation operators respectively to the right and the left of the T-product getting 



d 3 k 



(27r) 3 2w fc 



n(w k )(0\a(p', f)a"(k, f)T{Ci(xi) ■ ■ ■ h'fa) ■ ■ ■ £/"(^) • • • £/(zw)« + ft. f)<J(p, /)|0) paircd , (3.10) 



where the subscript n denotes the n-dependent contribution from the ^-propagator connecting Xi and Xj . 

It remains to show that we do not get any additional sign, if the x field is fermionic. First note that bringing 
the fields x( x i) ancl X.( x j) together to form the propagator and then putting them back to their old positions can be 
effected by the same set of field interchanges , one in reverse order of the other. Thus we do not encounter any extra 
minus sign here. Also the interaction Lagrangians being quadratic in fermionic x-> the operators a(k, /') and <v(k, /') 
do not produce any minus sign while moving through them. However, in the left hand side of Eq (3.8) an initial 
interchange of a and is needed which produces a minus sign to cancel the minus sign in front of the first term in 
Eq.(3.6). Thus Eq.(3.10) remains valid for both fermionic and bosonic operators. 

Eq.(3.10) gives the linear contribution from a definite x-propagator, namely the one between Xi and Xj. To get the 

total contribution from all the x~ propagators in the diagrams, we must allow a(k, /') and (r(k, /') to be paired with 
X fields at all vertices in all possible ways, not just with x( x i) an d x( x j)- Thus we get 

/d 3 k -> 
^—^—n(Lu q )(0\a(p', f)a(k, f')T{d(xi) ■ ■ ■ dix^ja^k, f)a\p, /)|0) pairod . (3.11) 

This matrix element is just the sum of all Feynman diagrams in coordinate space in the N-ih order of perturbation 
expansion for the scattering amplitude T//< introduced previously by Eq.(2.7). We thus prove Eq.(2.8) in an arbitrary 
order of pertubation theory. 




FIG. 1. The virial expansion of self-energy in medium to first order. 
IV. DISCUSSION 

It is well-known that the effective theory incorporating the symmetries of the QCD Lagrangian, called the chiral 
perturbation theory, can successfully describe the strong interaction processes in the low energy region. This theory 



4 If /' is a fermionic species, there may be two such propagators between Xi and Xj. We, of course, have to extract this 
contribution from both of them. 
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finds a natural application in the realm of hadronic statistical physics [|3|-|T5[. There is also a phenomenological 
approach in statistical mechanics for interacting systems, namely the method of virial expansion. 

In a region where the expansions of both the methods converge rapidly, the virial expansion would prove to be an 
identity in chiral perturbation theory. However, there are situations, due to the proximity of resonances, for example, 
where the effective coupling constants in the chiral Lagrangian can be rather large and the leading term in the chiral 
expansion may hold only in a limited region of interest. The virial expansion, on the other hand, may enjoy a wider 
range of validity. 

The case of the nucleon self-energy at finite temperature illustrates this point Q. Due to the presence of A(1232) 
resonance near the 7T./V threshold, the chiral expansion converges slowly, whereas the first term in the virial expansion 
is a good representation for low enough pion densities. Another example is the nucleon self-energy in nuclear medium, 
which involves the interaction of the two-nucleon system. Here the presence of bound or virtual two-nucleon states 
close to the threshold of NN scattering makes it difficult to formulate a satisfactory chiral perturbation theory for 
this system |l6|] . But the virial formula is expected to hold for densities close to the nuclear saturation density |L7j . 

The proofs for the virial expansion described here are indeed simple, mainly because we restrict to the first order 
formula. But the methods are not restricted to the first term in any way. They may well provide simple alternative 
proofs in the real time formalism for the complete virial expansion Q. One has only to take care of two additional 
aspects, namely, the disconnected parts that result from opening more than one internal lines and the 2 ® 2 matrix 
structure due to the so-called ghost vertices. 
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